Heavy Quarkonium Production through the Semi-Exclusive e~^e 
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Within the framework of tlie non-relativistic QCD, we make a detailed discussion on the heavy 
quarkonium production at a e^e~ collider with high luminosity up to £ cx 10^*~^®cm~^s~^ and with 
the collision energy around the peak. They are produced through the semi-exclusive channels 
e"'"e~ — >■ \Hqq) + X with X = QQ or gg. Here the heavy quark Q = b or c respectively. It is 
noted that in addition to the color-singlet IS-level quarkonium states, the color-singlet quarkonium 
states at the 2S'-level and the IP-level, and the color-octet quarkonium state \{QQ)[^^ Sf'^g) can 
also provide sizable contributions. The heavy quarkonium transverse momentum and rapidity 
distributions for the e+e~ collision energy Ecm ~ rnz are presented. For the charmonium or 
bottomonium production via the Z'^ propagator, there is approximate "spin degeneracy" for the 
spin-triplet and spin-singlet quarkonium states. Furthermore, uncertainties for the total cross 
sections are predicted by taking = 1.50 ± 0.15 GeV and rUb — 4.90 ± 0.15 GeV. Around 
Ecm = rnz, due to the Z°-boson resonance effect, total cross sections for the channels via the 
Z°-propagator become much larger than the channels via the photon propagator. Then, in addition 
to the B factories as BABAR and Belle and the hadronic colliders as Tevatron and LHC, such a 
super ^-factory will provide an excellent platform for studying the heavy quarkonium properties. 



PACS numbers: 13.66.Bc, 12.38. Bx, 12.39. Jh 



I. INTRODUCTION 

The investigation of the heavy quarkonium, e.g. char- 
monium and bottomonium, is helpful for us to achieve 
a deeper understanding of QCD in both the perturba- 
tive and nonperturbative sectors A comprehen- 

sive review of heavy quarkonium physics can be found 
in Ref.[l2j- In comparison to the hadronic colliders as 
Tevatron and LHC, an e"'"e~ -collider has many advan- 
tages, such as it provides a cleaner environment and the 
collision energy of the incoming electron-positron beams 
is well under control. Recently, a super Z factory run- 
ning at an energy around the Z^-boson mass with a high 
luminosity C ~ 10'^*~'^®cm~^s~^, similar to the GigaZ 
program of the Internal Linear Collider [l^, , has been 
suggested In the present paper, we will concentrate 
our attention on the heavy quarkonium production at 
such a super Z factory. This can be a useful reference for 
experimental studies and also be a useful compensation 
for the present BABAR and Belle results on the heavy 
quarkoniums. 

Within the framework of non-relativistic quantum 
chromodynamics (NRQCD) [l6|, the production process 
can be factorized into a sum of products of the short- 
distance coefficients and the NRQCD long-distance ma- 
trix elements. The short-distance coefficients are per- 
turbatively calculable in a power series of as at the 
energy scale around the heavy quark mass, which can 
be determined by matching the results obtained under 
full QCD and NRQCD accordingly. Generally, the non- 



perturbative long-distance matrix elements can be de- 
termined from lattice QCD calculations, or by fitting the 
prediction with the experimental data, or be roughly esti- 
mated by means of NRQCD velocity power counting rule. 
It has been shown that the matching of the color-octet 
matrix elements from the hadronic experiments strongly 

tends on the parton distribution function (PDF) (iTI - 
and the PDF uncertainty usually provides one of the 
key errors for theoretical estimation. Regarding to this 
point, the e"*"e~ collider provides a better platform for 
precise studies than the hadronic colliders. 

Within the framework of NRQCD, the physical state of 
a heavy quarkonium is generally regarded as a superpo- 
sition of Fock states, and the relative importance among 
those infinite ingredients is evaluated by the velocity scal- 
ing rule [l6l]. When the {QQypaii in the hadron l-f^gg) 
is in color-singlet state, it usually gives the dominant 
contribution to the heavy quarkonium productions and 
decays; while the color-octet {QQ)-pa,u: is usually powers 
of V suppressed. Here v stands for the typical velocity of 
the heavy quark and anti-quark in the quarkoniom rest 
frame, v'^ ~ 0.23 for (cc)-quarkonium and ~ 0.08 for 
(66)-quarkonium. For example, in the velocity expansion, 
we have 



\hQ) = 0(1)|QQ[^P(^)]) 
\XQj) = O(l)IQQpPj')]) 



-Oiv)\QQ['S^,%) 
-Oiv)\QQ['s[%) 
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where Q — b or c. Throughout the paper, we denote 
the ((5Q)-quarkonium with an extra superscript (8) to 
stand for the color-octet state and the superscript (1) 
to stand for the color-singlet state. Later on, without 
introducing any confusion, we omit the superscript (1) 
for the color-singlet state. The angular momentum of 
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the quarkonium is set in the square brackets. The color- 
octet {QQ)-pa.ii can give sizable and observable contribu- 
tions in certain cases or in certain kinematic regions when 
the color-singlet terms are highly suppressed by the hard 
scattering part. For example, the energy spectrum of the 
charged lepton for the color-octet components of the Be 
meson is dominant over its color-singlet component when 
the lepton has high energy [22,] . 

In the present paper, we will make a detailed discus- 
sion on the heavy quarkonium production at the super 
Z factory via the following two semi-exclusive channels: 

I 



and 



where the heavy quark Q — c or b respectively. It is 
noted that due to the absence of Fig.(j4h) (by Furry The- 
orem) for the ipQ production from the channel e~^e~ — 
— > -i/jQ + QQ, the contribution of the component 

l(QQ)[(l'^-Pj^^)]ff) is negligible. In this paper, we will 
only focus on the dominant color-octet channels listed 
above ^. 

Principally, there are two approaches to deal with the 
heavy meson hadroproduction. One is the fragmentation 
approach, which automatically sums up the dominant 
contributions, including some important higher-order ef- 
fects, into the total/differential cross sections by using 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution 
equation. The fragmentation approach is comparatively 



^ There are other less important channels, for example, for the 
channel e+e~ — > \Hqq) + gg, we find its total cross-section from 
the color-singlet P-wave quarkonium states is also negligible. We 
will not consider them in the present paper. 



e+e- 7*/Z" \Hqq} + QQ and e+e" ^ 7*/^° ^ 
\HQQ)+gg. At the present, these two processes, together 
with the present BABAR and Belle data, are widely 
adopted for estimating the color-octet components. How- 
ever, the present estimations of them (especially for the 
charmonium case) are inconsistent with each other, c.f. 
Refs.[23l-[25j. Thus, it would be helpful to find a new 
platform, such as the super Z factory, to learn more of 
its properties. More explicitly, we will deal with the fol- 
lowing production processes: 



simple and can reach to leading logarithm order or even 
higher orders more easily, but is satisfied only in the cases 
if one is only interested in the produced meson itself. The 
other one is the so-called complete calculation approach, 
in which we directly deal with the full hard scattering 
amplitude without any approximations. This approach 
sometimes is hard to derive the analytical expression even 
at the leading order, but this approach retains the in- 
formation of the accompanied quark or gluon jets and is 
experimentally more useful. In the present paper, we will 
mainly adopt the complete pQCD calculation approach 
to deal with these channels. 

In comparison to the B factories as BABAR and Belle, 
we will show that a large number of heavy quarkonium 
events can be generated due to the Z°-boson resonance 
effect. Some interesting examples for the heavy quarko- 
nium production at such a super Z factor can be found 
in Refs. (26l - l32| . Moreover, it is found that at the super 
Z factory, some novel features will appear. For exam- 
ple, when (QQ) is in color-singlet, the channel e~'"e" 
Z^ — >■ \Hqq) -f QQ can not distinguish the spin-singlet 



e+e- ^ Z°,7* ^ |(QQ)[(l'5i) , {T" S^) , (l^So) , (21^0)]) + QQ, 
e+e- ^ ^0,7* ^ |(QQ)[(liPi) , (l^Po) , (l^Pi) , (l^Pa)]) + QQ, 

e+e- ^ ^ |(QQ)[(li4^') , (l^^f )].g) + QQ ^ i^q + QQ 
e+e- ^ ^ |(QQ)[(li5^^^)]5) + QQ ^ Hq + QQ 
e+e- ^ ^ \{QQ)[(l's[^^)]g) + QQ ^ XQJ + QQ 

e+e- ^ Z°,r ^ \{QQ)[{l^Si) , {2^3^)]) + gg, 

e+e- ^ ^ |(QQ)[(l'^o) , {2'So)])+gg, 

e+e- -^Z'^ |(QQ)[(li5r) , (l=^5("^)].g) + g ^ i^Q + g 

e+e- -^Z°^ |(QQ)[(li4"))].g) + g ^ hg + g 

e+e- -^Z°^ \iQQ)[(l's[^^)]g) +g^ XQ.i + 9 
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and spin-triplet S'-wave states, which can be regarded as 
the approximate "spin degeneracy" . 

The remaining part of the paper is organized as fol- 
lows. In Sec. II, we present the calculation technique for 
dealing with the heavy quarkonium production processes 
e+e- \Hqq) + QQ and e+e" \Hqq) + gg, where 
the intermediate {QQ)-state is in either color-singlet or 
color-octet state respectively. In Sec. Ill, we present our 
numerical results and discussions. Total and differential 
cross sections are discussed, and we adopt the fragmen- 
tation approach for a further explanation of the "spin 
degeneracy" . Sec. IV is reserved for a summary. 



II. FORMULATION AND TECHNIQUE 



According to the NRQCD framework, the differential 



cross section for the process, e+e 
factorized as [H, HI] : 

da = do- (e+e- ^ {QQ)[n] 



QQ/ 



X , can be 



^\ (0|O^(«)|0) /o^ 



The production matrix element (010^(71)10) is propor- 
tional to the inclusive transition probability of the inter- 
mediate perturbative {QQ)-j>a,iT in [n]-state into the final 
bound-state \Hqq). The symbol [n] = [m^^+'^L^P'^^'>] 
with definite energy level m, spin S, orbital angular mo- 
mentum L and total angular momentum J, represents 
the above mentioned states for the intermediate (QQ)- 
pair respectively; i.e., 

n = l\So, 2' So, I'Pi, 1^5^^^; l^^i, 2^51, l^Pj, l''s[^\ 

with J — 0,1,2. These states provide the dominant con- 
tributions to the processes e+e~ 7*/Z° \Hqq) + 
QQ and e+e- ^ 7*/^° ^ \HQQ)+gg up to 0{v^). The 
parameters Ncoi and Npoi refer to the number of colors 
and polarization states of the intermediate (QQ)-pah:. 
Ncoi = 1 for the color-singlet state or Ncoi — — 1 
for the color-octet state. As for the color-singlet compo- 
nents, their matrix elements can be directly related to 
the wavefunction at the origin or the first derivative of 
the wavefunction at the origin, which can be computed 
via the potential modes and/or potential NRQCD and/or 
lattice QCD. As for the color-octet matrix elements, they 
can be estimated by using the velocity scaling rule or be 
determined experimentally. 

The short-distance cross-section da{e^{p2)e^ (pi) — > 
{QQ)[n] + X) can be written in the following form: 



{QQM 



X) = 



4\/(pi ■P2y - 



(4) 



where k stands for the number of final state particles, ^ 
means that we need to average over the spin states of 
the initial particles and to sum over the spin and color 



of all final particles. The phase space with k final-state 
particles : 



d$fc - (27r)*<5* \pi+P2 




k+2 



J-^i (2^)32// 



The phase-space integration can be done with the help 
of a combination of the subroutines RAMBOS [s^ and 
VEGAS [35[, which can be found in generators GENX- 
ICC [ll and BCVEGPY After generating proper 

phase-space points, one can calculate the numerical value 
for the squared amplitude |A^p. Furthermore, in deriv- 
ing the fragmentation functions, a splitting phase space 
is helpful. To be useful reference, we put some details for 
doing the phase-space splitting in the appendix. 

The hard scattering amplitude for those processes can 
be written as, 



(5) 



(6) 



The leptonic current 

L';.,,^Vrip2)r^Ur>ipi). 

The value of i is process dependent, e.g. 

1 = 4 for e+e- ^ {QQ)[n] + QQ 

i — 6 for e+e- {QQ)[n] + gg 

for {QQ)[n] in color-singlet state; while for the case of 
color-octet production, 

z = 2 for e+e- ^ \{QQ)[(l' S^,^^)]g) + g 

for e+e- ^ ^ \{QQ)[{l' s[^^)]g) + g 

I = 6 for e+e- ^ ^ \{QQ)[(l^ sl^^^]g) + QQ 

i = 8 for e+e- ^ ^ |(OQ)[f I'^^'lls) + QQ 



The spin-indices r and r' stand for the spin projections 
of the initial electron and positron. For the quarkonium 
production through the Z'^-boson propagator, the vertex 
— 7^(1 — 4sin^ 9yj — 7^) and the propagator Z?^i, = 
k^-ml+imzr, i'^d^iiy + kf.k^/k'^) , where F^ stands for the 
total decay width of Z° boson. While, for the quarko- 
nium production through the photon propagator: the 
vertex F^ = 7^ and the propagator I?^^ = —-^g^v 

The overall parameter C — Cs or Co for the production 
of the color-singlet or the color-octet quarkonium state 
respectively. The simplified hard scattering amplitudes 
A'( are process dependent, which will be presented in the 
following subsections. 



A. Color-Singlet Case 

1. e+(p2)e-(pi) ^ \HQQ){pz) + Q{pi)Q{ps) 

Typical Feynman diagrams for the process e+e~ — >■ 
\Hqq) -\- QQ through color-singlet (QQ)-quarkoniimr 
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.QiPr,) 




Qipi) 



Q{pi) 



+ 2 flipped oraphs 



states are presented in Fig.([T]). According to the Feyn- 
man diagrams, we can write down the simphfied hard 
scattering amphtudes of the short-distance cross- 
section (i(7(e+e~ — ^ {QQ)[n] + QQ). 



FIG. 1: Feynman diagrams for e'^e~ — > \Hqq) + QQ via the 
perturbative state (QQ)[n\- Here Q stands for c or b quark, 
and [n] stands for the color-singlet quarkonium state. 



For \n] stands for the color-single S'-wave states, we 
have 



J 



A ="s(P4)Fqq- ^ 2-7p- y^lpVs'(p^), 

(P3 + PbY - (p3i + p^r 
(Pi + P5 + P3ir - ^^[P3i+P5r 



{P32 + PiY {P3 + PiY ~ m^Q ' 
ttO(1) / \ 

^W^p. -fi-h-f32^'^Q ^ 

(P32 + P4)2 ' (P4 + P5 + P32)^ - m| ^"^^ ^^^^ 



(7) 
(8) 

(9) 
(10) 



r 



For stands for the color-singlet P-wave states, we have 
I 



and 



Ai ^ Us[pi)e^{p:i) — 



'^'^{P3+P,?-ml^'{p3i+P,)'' 



■'(Pb), 



q=0 



ji',S=0,L=l 



^i',S=0,L=l 



= usiPiya(.P3)-r- 

dqa 



7p 



.(P4).L(P3)^ 



A^ =Us{pi)eJyP3)-7— 



(p4+P5+P3l)2-m| QQ(;,3i+p5)2> 

^'{p32+P4r^'{p3+PA?-ml W 

-n-n \ , , , . 



9=0 
Vs'iPd), 



q=0 



n(QQ)fe) 

7p 



iWr-- p. ^I^4-|^5~^32 + "^Q 
(P32+P4)2 '5'5(p4+P5+P32)2-m|^'' 



9=0 



i/,S=l,L=l - ^„ ^_J,- / s d 



|^3-|^5+"1Q n^QQ)b3) 



pf » 3 



{P3 + PbY - m'^ {p3i + ps) 



Tip 



Vs' {P5 



q=0 



A^ 



v,S=l,L=l 



A 



i/,S=l,L=l 



UsiP'i)£al3iP3) 



uAPi)ei'piP3) 



dqa 

d 

dqa 



4 ' ^5 ' ^31 



(P4 + P5 + P3l)^ - m2 ' W (p3i + p5)2 



^''(p32+P4)2^''b3+P4)2-m^ W 



9=0 



9=0 



(11) 
(12) 

(13) 
(14) 



(15) 
(16) 
(17) 
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A. 



u,S=l,L=l 



n 



(QQ) 



"32 



{P32 + PiY (P4 +P5+ P32Y - ml 



-Ip 



Vs'{p5) 



(18) 



9=0 



The parameters s and s' stand for the spin projections 
of the outgoing quark and antiquark respectively. 

For convenience, we have introduced a general interac- 
tion vertex 



where Jz = Sz or Iz respectively. For the spin-triplet 
P-wave states {^Pj with J = 0,1,2), the sum over po- 
larization is given by (ssj 



(19) 



(1 - 7'^)/2 and Pr = {l+j'^)/2. Here a 
and (3 = —^s'm^Ow for (Zcc)-vertex, a 



where Pl 

4 

and /3 = 1 for (7*(50)-vertex, respectively. 
The momenta of the constituent quarks are 



fsin^ 

2—4 sin^ On, and /3 = — | sin'^ 9w for (Z66)-vertex, a 



P31 



Mqq 



P3 + q and 



Mqq 



P3 - q, 



(20) 



where Mqq = niQ + ttiq is implicitly adopted to ensure 
the gauge invariance of the hard scattering amplitude, q 
is the relative momentum between the two constituent 
quarks inside the quarkonium. 

The covariant form of the projectors are 



QQ , 



■irnqmr 



'^32 



"^q)75(|>3i +™q) (21) 



and 



with 



(22) 



n(QQ)(P3) 



ArnqniQ ^2 



Inserting these projectors into the amplitude, the ampli- 
tude can be squared, summed over the spin in the final 
state and averaged over the ones in the initial state. The 
selection of the proper angular momentum is done by 
performing proper polarization sum. We define: 



(24) 



The sum over polarization for a spin-triplet S'-wave state 
i'^Si) or a spin-singlet P-wave state (^Pi) is given by: 



(25) 




5 permutations 



e+(P2) 



FIG. 2: Typical Feynman diagrams for e^e~ — >■ \Hqq) + gg 
through the perturbative state {QQ)[n], here Q stands for the 
c or 6 quark and [n] represents the color-singlet quarkonium 
states. 



Jz 



Jz 



For the overall color-singlet parameter, we have Cs — 
iTf sin^ eJ(tcose^)^ ^ij quarkonium production 

through Z*'-boson propagator and Cs = '^/^^Q^'^ 9s^ij ^^r 
the quarkonium production through photon propagator, 
where eg stands for the electric charge of Q in unit e. 



2. e+{p2)e (pi) ^ \Hqq){p3)+ g{p4)gip5 



Typical Feynman diagrams for the process 



\H. 



QQ/ 



gg are presented in Fig. ([2]). According to the 
Feynman diagrams, one can write down the simplified 
hard scattering amplitudes A'j^ for the short-distance 
cross-section d(T{e~^e~ — {QQ)[n] + gg) as 



Tr 



KIM, 



QQ 



'^32 



mq 



{P32 +P4+P5 y - ml 



iip^) 



'^32 



-f^ + mq 



{P32 +PiY 



/(P4) 



(26) 
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and the amplitudes A2 , ■ • • , Aq can be derived by per- 
mutations. 

For the overah color-singlet parameter, we have Cs = 
i73 sin^ eJitcose^y ^<^b for the quarkonium production 
through Z°-boson propagator and Cs = ^^eqe^g^Sab 
for the quarkonium production through photon propaga- 
tor, where a,b — (1, • • • , 8) are color indices of the two 
outgoing gluons. 

The simplified amplitudes for the case of color-singlet 
P-wave states can be written down in a similar way as the 
process e~^e~ — >■ \Hqq) + QQ. Numerically, it is found 
that the total cross-sections for the color-singlet P-wave 
states are negligible, so we do not present them here. 

B. Color-Octet Case 

1. e+{p2)e-{pi) ^ Z° ^ 

, (l'S(^')]5)(P3) + Q(P4)Q(P5) 

Typical Feynman diagrams for the processes e+e^ — > 
Z° ^ \iQQ)[{l^S^^^^]g) + QQ and e+e" ^ Z° ^ 



^5 = Ms(p4)7pTr 



• The remaining two simplified amplitudes A^ and 
A^ for the case of \{QQ)[(l^s[^^)]g) production. 



For the overall color-octet parameters, we have Co — 



\{QQ)[{l^Sr'j]g} + QQ are presented in Figs.(I311) ac- 
cordingly. It is found that 



• Four Feynman diagrams as shown in both Fig. ([3^) 
and Fig.(|3^) have the same topologies as those of 
Fig-©- Their simplified amplitudes A^, - - - , A4 are 
the same as those of Eqs. (|7l8l9ll0p . While, the 
overall parameter Cs there should be replaced by 
the present case of Co, and the color-singlet ma- 
trix element should be replaced by the color-octet 
matrix element. 



• The two Feynman diagrams as shown in Fig.(j3l3) 
and Fig.(l4}3) have the same topologies, and their 
simplified amplitudes A'^ and Aq are the same. 
More explicitly, we have: 



Vs'{P5), (27) 



Vs'{p5). (28) 

I 

as shown in Fig. ([4};), are 



(29) 
(30) 

I 

color factors cf"^'^'''''''^'', according to Figs.(j4^|4j3|4t), are: 



= Y.(T%.- X iV2T%„, X (T'^ymj (31) 

m,n 



{Pi+Pby (P32 +P4 +P5)^ - "^Q 

{P31 + P4 + Ps)'^ - m% ^^{pi+p^f 



{Pz + P^f - m^Q 

"njQQ)(P3) 



n?oo^(P3 



2 IP 
PI 



Vs'iPz), 



"6(p4)7pTr 



2 -7p 
Pi 



2— ;;rrQQ«s'(P5)- 



(P3 +PiY - mi 
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_yQ[l'5,n(p;,)g-(p^j 




.i„f81, , , 



- 2 flipped graphs 



■ Qipr.) + 1 flipped graph 



FIG. 3: Typical Feynman diagrams for e^e — > - 
|(QQ)[^1^5Q^''^]g) + QQ, where Q stands for c or b quark. 



e (Pi) 



oO[(i'5,rj,(i-'sS'*')](P3 




+ 1 i)cniiutatioii 



FIG. 5: Typical Feynman diagrams for e'^e — >■ Z° — >■ 
l(QQ)[(l^S*^') , (l^5'J^^)]5> + 5, where Q stands for the 
heavy c or 6 quark. 



(a) C'^^'Z' 




(?0[r'Sr'](p3) 

+ 2 flipped graphs 



■ Qipi] 

■ 0(P5) 

i 00[i'sr'l(p3) 

Q(P5) 

I eO|i''s!"'](p,) + 1 flipped graph 

FIG. 4: Typical Feynman diagrams for e^e~ — )■ — >■ 
|(QQ)[(l^5f + QQ, where Q stands for c or fe quark. 

= ^ (^")-" ^ (^7^'')n™ (32) 

= ^ (^")"" ^ (^^r'^)"™ (33) 

where i, j, m,n — (1, 2, 3) are color indices of the outgoing 
anti-quark Q, quark Q and the two constituent quarks 
Q and Q in the heavy quarkonium, respectively. The 
superscript indices (a), (&) and (c) of refer to the 
corresponding figures in Figs. pi4p . After simplification, 
we obtain 

(a) ^ _ V2.T.d ^(b) ^ ^(c) ^ 

Here d — {1, ■ ■ ■ ,8) represents the color of the color-octet 
quarkonium. 

2. e+(p2)e-(pi) ^ Z« ^ 
|(0Q)[(l^5(«'),(l^5f))l3)(p3) + 3(P4) 

Typical Feynman diagrams for the processes e+e^ — > 
^ |(QQ)[(liS'^^^) , (l^S-^^^^Jg) + g are presented in 
Fig.®. Its first simplified amplitude takes the form 



^5' = Tr 



'^32 



and the second amplitude A2 can be obtained from Ai 
by permutation. 

For the overall color-octet parameter, we have Co ~ 
^ X bad- Here a,d = (1,- • 



2 sin^ e„(4cose„)2 



represent 



the color indices of the outgoing gluon and the color octet 
quarkonium, respectively. 



III. NUMERICAL RESULTS 



Matrix element 


|*(0)|"or|*'(0)|^ 




1.2 GeV'* 


0.067 GeV=* 


(OlOf (^5i)|0) 


7.6 X 10"' GeV^ 


0.042 GeV=* 


(0|C?(3Si)|0) 


9.3 GeV^* 


0.517 GeV^* 


(oior('si)io) 


4.6 GeV'* 


0.256 GeV=* 


{0|O?=i(3Pi)|0) 


3.2 X 10"^ GeV'' 


0.018 GeV^ 


{0|O^" (3Pi)|0) 


6.1 GeV^ 


0.339 GeV^ 


(0|Og/'^(='Si)lO) 


1.2 X 10"^ GeV^ 




(0|O8^(»Sl)|0) 


9.5 X 10"^ GeV^ 




{0|O3^-(«5i)|0) 


1.6 X 10-^ GeV^ 




(0|O3^"(^5i)|0) 


4.3 X 10"^ GeV^ 





TABLE I: Input color-singlet and color-octet matrix ele- 
ments [3^, where the color-singlet matrix elements are related 
to the wavefunctions at the origin or the first derivative of the 
wavefunctions at the origin respectively. 

In doing numerical calculation, the heavy quark masses 



are taken as: = l-50+^:i^ GeV and mt = 4.90l[J i^ 
GeV. The color-singlet and color-octet non-perturbative 
matrix elements are taken from Ref . (38| . which are pre- 
sented in Table m As for the production matrix elements, 
we adopt the relation 

(0|0^,«je5i)|0)^3(0|0^,,](i5o)|0). 

The production color-singiet matrix elements and the 
wavefunctions at the origin are related by [l6l [33| : 



(0|O^=|0) (0|Oi/'^|0) 



(34) 



6iV, 

(0|gf |0) 
6A^c 



- 1*15(0)^ 
|*2S(0)P 



(35) 
(36) 
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and 



Of |0) (0|OfM0) (0|Of^|0) 



(37) 

As for the production color-octet matrix elements, we 
adopt the relation [16 1 



(0|Or('^i)|0) = ^(0|0^(i5o)|0) 



(0|Or('5i)|0) = ^(0|0^(i5o)|0) 



(0|Or('^i)|0) = ^(0|0^(i5o)|0) (38) 



Here the heavy-quark spin symmetry [16| has been im- 
plicitly adopted. That is, we will not distinguish the 
wavefunctions at the origin for the spin-singlet (^Sq) state 
and the spin-triplet (^6*1) state at the same energy level. 
For the bottomonium production color-singlet and color- 
octet matrix elements they have the similar relations as 
those of the charmonium. 

As for the renormalization scale, we take it to be 
2mc for charmonium production, and 2mb for bottomo- 
nium production. Other input parameters are taken 
from the Particle Data Group 0]: Tz = 2.4952 GeV, 
mz = 91.1876 GeV, sin26l^=0.2312. By using leading- 
order as running and as{mz) = 0.1184, we obtain 
asi2mc) = 0.237 and a^(2mf,) = 0.175. 



A. Properties of the color-singlet processes 

The color-singlet components provide the dominant 
contributions to the heavy quarkonium production pro- 
cesses. We will make a detailed discussion on the prop- 
erties of the color-singlet processes, which includes their 
total cross sections, their differential cross sections, and 
the uncertainties caused by the variation of the e~^e~- 
collision energy being slightly different from the Z^-boson 
mass, by the variation of the heavy quark masses and etc. 



1. Total and differential cross sections 

Total cross sections for the charmonium and bottomo- 
nium productions are presented in Tables |lT] and IIIII 
By adding all the P-wave charmonium or bottomonium 
states together, we obtain 

cr(e+e~ ^ \Hac) (IP) + cc) ~ 1.08 pb 



and 



\Hi,i){lP) + bb) ~ 3.25 X 10~2 pb. 



Furthermore, we have |5'25'(0)P/|^is(0)p ~ 0.6 for char- 
monium and ~ 0.5 for bottomonium. Then, in addition 
to the IS'-level quarkonium states, both the 2S'-level and 





Total cross sections 


production channels 


.Z^-propagator 


7*-propagator 


e+e" — i> ric{lS) + cc 


1.76 


4.09 X 10"^ 


e+e" ^ ri'rX2S)+cc 


1.11 


2.59 X 10"^ 


e+e" ^ J/ipilS) +CC 


1.83 


4.39 X 10"^ 


e+e- ■,p'{2S) + cc 


1.16 


2.78 X 10"^ 


e+e- ^ J/i>{lS) + gg 


3.84 X 10"^ 


6.99 X 10"'' 


e+e" V'(2S) + .g.g 


2.43 X 10"^ 


4.42 X 10"" 


e+e" rjc{lS)+gg 


1.72 X 10"^ 





c+e- ^n'c{'2S)+gg 


1.09 X 10"^ 





e+e" hc(lP) + cc 


2.34 X 10"^ 


5.60 X 10"'' 


e+e" -s- Xco(l-P) + cc 


3.32 X 10"^ 


7.78 X 10"" 


e+e" Xci(l-P) + cc 


3.66 X 10"^ 


8.37 X 10"" 


e+e" -5> Xc2{lP) + cc 


1.44 X 10"^ 


3.24 X 10"" 



TABLE II: Total cross section (in unit: pb) for the color- 
singlet charmnoium production at the ^"-peak [\fS = mz)- 
The channels are through and 7* propagators respectively. 
mc=1.5 GeV. 





Total cross sections 


production channels 


-propagator 


7*-propagator 


e+e" ribilS) + bb 


1.94 X 10"^ 


8.69 X 10"^ 


e+e- r]'b(2S) + bb 


9.60 X 10"^ 


4.30 X 10"^ 


e+e- -> T(IS') + 66 


2.18 X 10"^ 


1.16 X 10"" 


e+e" ^ r'{2S) + bb 


1.08 X 10"^ 


5.74 X 10"^ 


e+e- ^rilS)+gg 


6.54 X 10"^ 


9.16 X 10"^ 


e+e- -^T'{2S) + gg 


3.24 X 10"^ 


4.53 X 10"^ 


e+e- 776(15) +55 


8.18 X 10"^ 





e+e- ^4{2S)+gg 


4.05 X 10"^ 





e+e- ht{lP) + bb 


6.84 X 10"^ 


3.36 X 10"*^ 


e+e" ~^Xbo{lP) + bb 


1.09 X 10"^ 


5.18 X 10"*^ 


e+e- -> Xbi{lP) + bb 


1.04 X 10"^ 


4.22 X 10"*^ 


e+e" ^Xb2{lP) + bb 


4.30 X 10"^ 


1.67 X 10"*^ 



TABLE III: Total cross section (in unit: pb) for the color- 
singlet bottomonium production at the Z^-peak {s/S = mz)- 
The channels are through Z" or 7* propagator respectively. 
mi,=4.9 GeV. 



IP- level quarkonium states can provide sizable contribu- 
tions to the production channel e+e" — > \Hqq) + QQ. 

More explicitly, to show the relative importance of 
those channels via the Z'^ propagator, we define two type 
of ratios 



and 



cr„+ 



(39) 



(40) 



We obtain, 

^Vc+gg = 9.4% , Rr^'^+gg — 6.0% , Rj/^+gg — 2.1% , 
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Tl +CC 

'c 

J/\|;+cc 





E (GeV) 

cm * ' 



40 60 
E „ (GeV) 



100 



FIG. 6: Total cross sections of the channel e+ + e 
versus the e^e~ collision energy Ecm ~ ^/S. 



\Hcc) + X, including the S wave and P wave charmonium states, 




20 40 60 80 
E „ (GeV) 



100 




20 40 60 80 100 
E (GeV) 



FIG. 7: Total cross sections of the channel + e — > — > \Hcc) + X, including the 5* wave and P wave charmonium states, 
versus the e^e~ collision energy Ecm ~ The two curves for tjc + cc and J/ip + cc are almost coincide with each other, and 
the two curves for XcO + cc and Xci + cc are almost coincide with each other. 



Tl^+bb 
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cm * ' 
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E.„ (GeV) 
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FIG. 8: Total cross sections of the channel + e 
versus the e^e~ collision energy Ecm ~ ^/S. 



\Hf£} + X, including the 5* wave and P wave bottomonium states. 
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cm ^ ' 




100 
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FIG. 9: Total cross sections of the channel e"*" + e — )■ — +X, including the S wave and P wave bottomonium states, 
versus the e^e~ collision energy Ecm = y/S- The two curves for r]b + bb and T + bb are very close to each other, and the two 
curves for Xbo + bb and xti + bb are very close to each other. 



R-Xal+CC 

and 



1.3% , 

63% , Rh^^ 



: 96% , R 
13% , R 



XcQ-\-CC 



61% , 

= 18% 



20% , Rxc2+cc — 7.i 



R 



Vb+gg 



38% , Rjj'^+gg — 19% , Ry+gg ~ 30% 



R.T- 



bb 



R 



Xbi+bb 



15% 

50% , ii-hh+bb 

4.8% , R 



Rvb+bb - 89% , Rni+bb — 44% , 

R 

•"b I "-"^ 

bb = 2.0% 



- 3.1% , R^^g+hb - 5-0% 



Xb2- 



It has been observed that the Z°-boson resonance ef- 
fect will raise the production rate by several orders of 
magnitude in comparison to the Belle and BABAR ex- 
periments, then it could be a better platform for studying 
the properties of the heavy quarkoniums. 

For the e"'"e^ collision energy Ecm = = mz, the 
cross sections for the channels via the 7* propagator are 
much smaller than those channels via the propagator. 
For example, we have 



= 2.4 X 10"^ 



and 



<^e+e 



>Z'>^J/ip+cc 



*7 



tl±^ = 5.0 X 10-4. 



^Z°^r+bb 



To show this point clearly, we present the total cross 
sections for the color-singlet channels versus the collision 
energy Ecm in Figs. (j6l7l8|9p . In the low collision energy 
region, the production cross sections are dominated by 
the cases via the 7* propagator: main contributions are 
around 6 — 20 GeV for charmonium production and 10 — 
40 GeV for bottomonium production. While, for the case 
of propagator, it is small in low energy region, but 
has a peak value at Ecm = fnz due to the Z^-boson 



resonance effect. This shows that the Super-Z factory will 
provide another platform for heavy quarkonium physics 
in addition to Belle, BarBar, CLEO and LHC colliders. 

For the charmonium production channels, from Fig. ([7]) 
and Fig.®, it is found that the total cross sections for 
e^e'^ — > — > J/'0 -f cc and e^e~ Z^ ^ rjc + cc, and 
the total cross sections for e+e" — > — > Xca + cc and 
e+e" — >■ — >■ Xci + cc are almost coincident with each 
other, especially for larger collision energy. The condi- 
tions for the bottomonium cases are similar. This shows 
there is approximate "spin degeneracy" for the produc- 



tion channels e+e Z 
when the collision energy Ec 



\H. 



QQl 



QQ. More clearly, 













J l^+cc 














0"e+e- 




-^rih+bb 






^r+bh 


O'e+e-- 




^Xbl+bb 




->Z0- 


^Xba+bb 



mz, we have 
96% 
91% 
89% 
95% 



(41) 
(42) 
(43) 
(44) 



We present the heavy quarkonium transverse momen- 
tum (pt)-distributions and rapidity (j/)-distributions for 
Ecm = mz in Figs. (|10lllll2ll3p . These figures show that 
either for the pj-distributions or for the y-distributions, 
the two curves for e~^e~ — >■ Z° — >■ J/i(j + cc and e+e" 
Z'^ rjc + cc, and the two curves for e+e" Z'-^ ^ 
Xco + cc and e'^e~ — >■ Z'' — >■ Xci + cc are also very close to 
each other. The cures for the bottomonium production 
are similar. 

As a cross check / explanation for the approximate 
"spin degeneracy" , we adopt the fragmentation ap- 
proach. Using the fragmentation approach the most im- 
portant/dominant higher-order effects can be included 
by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
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P, (GeV) p, (GeV) 

FIG. 10: The charmonium pt-distributions for the production processes, e^e~ — > — > \Hcc) + X, at the e^e~ colhsion energy 
Ecm = mz- The left is for X = cc and the right is for X = gg. 




P, (GeV) Pj (GeV) 

FIG. 11: The bottomonium pt-distributions for the production processes, e'^e~ Z" l^^bt) + X, at the e^e~ coUision 
energy Ecm = mz- The left is for X = bb and the right is for X = gg. 
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evolution equation, and the fragmentation approach gives 
a reasonable approximation to the full tree-level cal- 
culation as long as the transverse momentum of the 
produced heavy hadron is large enough 40i-l52|. Move 
over, it has been argued that the fragmentation approach 
may be extended to lower pT region under the FONLL 
scheme ,53.. ,54] . 



FIG. 14: Feynman diagrams for calculating the total frag- 
mentation probabilities for c J /il) and c ^ rjc.- 





Jim 

c(pi) 



In the large pt region, the probability of ^ J/tp or 
rjc is proportional and dominated by the fragmentation 
of c — > J/tjj or c ^ rjc- For the purpose, we need to 
deal with the Feynman diagrams shown in Fig. (|14[) . By 
taking the axial gauge, the amplitudes of Figs.(fT4bfT4H) 
are suppressed by a factor {rric/rnz) in comparison to 
Figs.([14]a[14l:), so they are neglected in deriving the frag- 
mentation probabilities for c — ?> J/ip and c ^ rjc- Fol- 
lowing the same procedure as shown by Ref. [sol ] , one can 
obtain the total fragmentation probabilities: 



<lzD^^j/^{z,imc) 



256 



a.(3rn.f X '^!f (Q)'^ r^-571n2 



/•i 256 



30 



l^is(0)P ^773 
-^^(— -371n2|, 



(45) 
(46) 



Because the ratio of the two fragmentation probabilities is 

dzDc->j/^(z,3m,)/ / AzDc^r^Xz, ^mc) ^ 97%. 



/ 

Jo 



/1189 „\ „ , /773 „\ „ This is consistent with Eq. dlT]) and clearly shows that 

' — 571n2 ~ 0.124 and 37 In 2 ~ 0.120, ^ — 



I 30 ) ' V 30 / ' ' there is spin- degeneracy for e+e at the high collision 
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Ptcut{GeY) 


5 


10 


15 




1.75(96%) 


1.55(85%) 


1.29(70%) 




1.11(96%) 


9.81 X 10-1(85%) 


8.17 X 10-1(70%) 


^J/ii{lS)+gg 


2.89 X 10"^ (75%) 


1.96 X 10-^(51%) 


1.30 X 10-^(34%) 


^i,'(2S)+gg 


1.83 X 10^^(75%) 


1.24 X 10-^(51%) 


8.23 X 10-^(34%) 


°"»7c(lS)+cc 


1.68(95%) 


1.46(83%) 


1.17(66%) 




1.06(95%) 


9.24 X 10-1(83%) 


7.41 X 10-1(66%) 


'^'7c(lS) + 99 


1.47 X 10"^ (85%) 


1.08 X 10-1(63%) 


7.51 X 10-^(44%) 


^v'r(2S)+gg 


9.31 X 10^^(85%) 


6.84 X 10-^(63%) 


4.75 X 10-^(44%) 


(^hc(lP) + cc 


2.22 X 10-^95%) 


1.96 X 10-1(84%) 


1.61 X 10-1(69%) 




3.17 X 10^^(95%) 


2.86 X 10-1(86%) 


2.42 X 10-1(73%) 


°"Xcl(lJ=) + cc 


3.46 X 10-1(94%) 


3.08 X 10-1(84%) 


2.56 X 10-1(70%) 




1.37 X 10-1(95%) 


1.24 X 10-1(86%) 


1.05 X 10-1(73%) 



TABLE IV: Cross section (in unit: pb) for the color-singlet charmonium production at the Z"-peak (\/5 — mz) with different 
Ptcut- The channels are through Z*^ propagator. mc=1.5 GeV. Here the percentages in the parentheses represent the ratios 
between the cross sections with various Ptcut and the one without Pt cut. 



Ptcut (GeV) 


5 


10 


15 


'^T{lS) + bb 


2.10 X 10-1(96%) 


1.89 X 10-1(87%) 


1.59 X 10-1(73%) 


<^r'{2S)+bb 


1.04 X 10-1(96%) 


9.35 X 10-2(87%) 


7.86 X 10-2(73%) 


0-T(lS) + 99 


5.90 X 10-2(90%) 


4.63 X 10-2(71%) 


3.37 X 10-2(51%) 


0"T'(2S) + 99 


2.92 X 10-2(90%) 


2.29 X 10-2(71%) 


1.67 X 10-2(51%) 




1.88 X 10-1(97%) 


1.68 X 10-1(87%) 


1.38 X 10-1(71%) 


<^7l'^^(2S)+bb 


9.30 X 10-2(97%) 


8.31 X 10-2(87%) 


6.83 X 10-2(71%) 


'^»;(,(lS) + 99 


7.85 X 10-2(96%) 


6.80 X 10-2(83%) 


5.35 X 10-2(65%) 


0"»,^(2S) + 99 


3.88 X 10-2(96%) 


3.36 X 10-2(83%) 


2.65 X 10-2(65%) 


'^/ih(lP) + bb 


6.60 X 10-^(96%) 


6.00 X 10-^(88%) 


5.00 X 10-^(73%) 


'^Xhn(lP) + 65 


1.05 X 10-2(96%) 


9.40 X 10-^(86%) 


7.85 X 10-^(72%) 


'^Xm(lP) + bb 


9.95 X 10-^(96%) 


8.95 X 10-^(86%) 


7.60 X 10-^(73%) 


'^Xm(1P) + 65 


4.16 X 10-^(97%) 


3.84 X 10-^(89%) 


3.33 X 10-=^ (77%) 



TABLE V: Cross section (in unit: pb) for the color-singlet bottomonium production at the .Z'^-peak (x/S = mz) with different 
Ptcut- The channels are through Z° propagator. mj,=4.9 GeV. Here the percentages in the parentheses represent the ratios 
between the cross sections with various Ptcut and the one without Pt cut. 



energy. Using these expressions, one can conveniently 
obtain the total fragmentation probabilities for — > T 
and — > ?7b by changing TOc into m;,. More over, the P- 



wave case can be derived through a similar way, whose 
detailed derivations can be found in Ref. [Ill , and we can 
obtain similar conclusions. 



In experiment, no detector can cover all the kinematics 
of the events, so only some of the He or Hf, events can 



be observed. That is, the events with small pt and /or a 
large rapidity y cannot be detected. So, the cross sections 
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ycut 


0.5 


1.0 


1.5 


f"j/i/.(lS)+cc 


7.00 X 10^^(38%) 


1.27(69%) 


1.60(87%) 


'^^'(2S) + cc 


4.43 X 10^^(38%) 


8.04 X 10-1(69%) 


1.01(87%) 


^J/ii{lS)+gg 


1.82 X 10"^ (47%) 


3.02 X 10-2(79%) 


3.58 X 10-2(93%) 


0"V'(2S)+S9 


1.15 X 10^^(47%) 


1.91 X 10-2(79%) 


2.27 X 10-2(93%) 


°"»7c(lS) + cc 


6.71 X 10^^(38%) 


1.23(70%) 


1.55(88%) 


(^vi.(2S) + cc 


4.25 X 10^^(38%) 


7.79 X 10-1(70%) 


9.81 X 10-1(88%) 


'^'7c(lS) + 99 


7.88 X 10'^ (46%) 


1.32 X 10-1(77%) 


1.58 X 10-1(92%) 


'^»7;.(2S) + 99 


4.99 X 10^^(46%) 


8.36 X 10-2(77%) 


1.00 X 10-1(92%) 


(^hc(lP) + cc 


8.97 X 10-2(38%) 


1.63 X 10-1(70%) 


2.05 X 10-1(88%) 




1.25 X 10^^(38%) 


2.28 X 10-1(68%) 


2.88 X 10-1(86%) 




1.39 X 10-1(38%) 


2.53 X 10-1(69%) 


3.20 X 10-1(87%) 




5.48 X 10-2(38%) 


9.97 X 10-2(69%) 


1.25 X 10-1(87%) 



TABLE VI: Cross section (in unit: pb) for the color-singlet charmonium production at the Z"-peak {\^ ~ mz) with different 
Vcut- The channels are through Z° propagator. mc=1.5 GeV. Here the percentages in the parentheses represent the ratios 
between the cross sections with various ycut and the one without rapidity cut. 



Vcut 


0.5 


1.0 


1.5 


'^T(lS) + bb 


9.40 X 10-2(43%) 


1.67 X 10-1(77%) 


2.06 X 10-1(94%) 




4.65 X 10-2(43%) 


8.26 X 10-2(77%) 


1.02 X 10-1(94%) 


0"T(lS) + 99 


3.54 X 10-2(54%) 


5.63 X 10-2(86%) 


6.38 X 10-2(98%) 


0"T'(2S) + 99 


1.75 X 10-2(54%) 


2.78 X 10-2(86%) 


3.16 X 10-2(98%) 




8.26 X 10-2(43%) 


1.49 X 10-1(77%) 


1.84 X 10-1(95%) 


^v'^(2S) + bb 


4.09 X 10-2(43%) 


7.37 X 10-2(77%) 


9.10 X 10-2(95%) 




4.23 X 10-2(52%) 


6.87 X 10-2(84%) 


7.92 X 10-2(97%) 


0"t,^(2S)+99 


2.09 X 10-2(52%) 


3.40 X 10-2(84%) 


3.92 X 10-2(97%) 


(^h,,{lP)+bb 


3.01 X 10-^(44%) 


5.29 X 10-^(77%) 


6.48 X 10-^(95%) 


'^Xboiin + bb 


4.60 X 10-^(42%) 


8.26 X 10-^(76%) 


1.02 X 10-2(94%) 


'^XM^-^P) + bb 


4.47 X 10-^(43%) 


7.93 X 10-^(76%) 


9.77 X 10-^(94%) 


^Xh■?.^-^P) + bb 


1.84 X 10-^(43%) 


3.27 X 10-^(76%) 


4.02 X 10-^(93%) 



TABLE VIL Cross section (in unit: pb) for the color-singlet bottomonium production at the .Z'^-peak (y/S = mz) with different 
ycut- The channels are through Z'^ propagator. m6=4.9 GeV. Here the percentages in the parentheses represent the ratios 
between the cross sections with various ycut and the one without rapidity cut. 



with different ptcut and i/cut are presented in Table IIVI 
TabltEl TablcEland TablfJVlIl 



2. Uncertainties caused by the variation of Ecm and mg 

For the leading-order calculation, the uncertainty 
sources include the bound-state matrix elements, the 
renormalization scale, the quark masses mf, and rric, and 
etc.. The conventional scale-setting assigns the typical 
momentum flow of the process as the renormalization 
scale, e.g. 2mc for charmonium production and 2mi, 
for bottomonium production. This ad hoc assignment of 
scale and its range leads to an important systematic error 
in the present theoretical estimations. In the literature, 
the principle of maximum conformality (PMC) (55| pro- 



vides a feasible way to derive a precise QCD predictions. 
The main idea of PMC is to sum all the non-confornial 
/3-terms in the perturbative expansion into the running 
coupling, and then the remaining terms are identical to 
that of a conformal theory. The PMC estimation is then 
scheme independent, and the remaining scale dependence 
can be greatly suppressed. At the present, the matrix 
elements and the strong coupling constant as emerge as 
overall factors and their uncertainties can be conveniently 
discussed when we know their values well, so we will not 
discuss their uncertainties in the present paper. 

As shown in the last subsection, if the collision energy 
Ecm is around mz, the total cross sections for the chan- 
nels via the 7* propagator are much smaller than those of 
the channels via the Z° propagator. In the present sub- 
section, we will adopt the channel via the Z°-propagator 
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Ecm 


97% mz 


mz 


103% mz 




6.66 X 10^^ 


3.84 X 10^2 


6.57 X 10^=^ 


<^^' + gg 


4.22 X 10~^ 


2.43 X 10"^ 


4.16 X 10"=* 




2.97 X 10"^ 


1.72 X 10"^ 


2.94 X 10"^ 


'^v'c+aa 


1.88 X 10"^ 


1.09 X 10"^ 


1.86 X 10"^ 




3.02 X 10"^ 


1.83 


3.26 X 10~^ 


^ip' -\-cc 


1.91 X 10"^ 


1.16 


2.06 X 10"^ 




2.91 X 10"^ 


1.76 


3.14 X 10"^ 




1.84 X 10"^ 


1.11 


1.99 X 10"^ 


ffh^ + cc 


3.88 X 10^^ 


2.34 X 10"^ 


4.18 X 10"^ 




5.51 X 10"^ 


3.32 X 10"^ 


5.92 X 10"^ 




6.08 X 10^^ 


3.66 X 10^^ 


6.55 X 10"^ 




2.40 X IQ-^ 


1.44 X 10"^ 


2.59 X 10"^ 



TABLE VIII: Total cross sections (in unit: pb) for the pro- 
duction channel e+e~ — > Z° — > \Hcc} + X. 





97% mz 


mz 


103% mz 


<^r + gg 


1.12 X 10"^ 


6.54 X 10~^ 


1.13 X 10-2 


^r' + gg 


5.54 X 10"^ 


3.24 X 10-2 


5.59 X 10-^ 


^Vb+aa 


1.40 X 10~^ 


8.18 X 10^2 


1.42 X 10-2 


'^v'^+aa 


6.93 X 10"^ 


4.05 X 10-2 


7.03 X 10-=^ 




3.57 X 10"^ 


2.18 X 10-^ 


3.95 X 10-2 


'^T'+bb 


1.77 X 10"^ 


1.08 X 10-^ 


1.95 X 10-2 


'^Vh + bb 


3.17 X 10"^ 


1.94 X 10-^ 


3.53 X 10-2 


'^ri'^ + bb 


1.57 X 10~^ 


9.60 X 10-2 


1.75 X 10-2 


<^H+bb 


1.13 X 10"^ 


6.84 X 10-^ 


1.24 X 10-^ 


'^Xhn + bb 


1.80 X 10"^ 


1.09 X 10-2 


1.95 X 10-=* 


'^XM+bb 


1.71 X 10"^ 


1.04 X 10-2 


1.88 X 10-=^ 


'^Xb-y + bb 


7.06 X 10^* 


4.30 X 10-^ 


7.74 X 10-'' 



TABLE IX: Total cross sections (in unit: pb) for the produc- 
tion channel e+g- — 5> Z° — > IH^i) + X. 



to do our following discussion, i.e. to deal with the chan- 
nel e+e" — T' — >• \Hqq) + X. For clarity, to discuss the 
uncertainty from any one parameter, we will take others 
to be their center values. 

To show the sensitivity of the total cross sections to 
the collision energy around the Z° peak, we calculate 
the total cross sections by taking Ecm = (1 i 3%)mz- 
Our results are presented in Tables I Villi and IIXI 

By varying Ecm within the range of (1 ± S%)mz, the 
total cross sections change greatly: e.g. it is down to 
(i7%) of its peak value for e+e" — )• Z° — )• J/^ + cc; while 
for the channel e+e- ^ T + bb, it is down to (11%)- 

Next, we discuss the uncertainties from the heavy 
quark masses by varying nic = 1.50 ± 0.15 GeV and 
nib — 4.9 ± 0.15 GeV respectively. Our results for 
e+e- Z° -J' \Hqq) + X are presented in Tables |X] 
and IXIl These tables show 

• For the charmonium production, when X = cc^ the 



TTic GeV 


1.35 GeV 


1.65 GeV 


uncertainty 


•^Jj^+aa 


4.55 X 10-2 


3.29 X 10-2 


+0.007 
-0.006 




2.88 X 10-2 


2.08 X 10-2 


+0.005 
-0.004 


^Vc+gg 


2.04 X IQ-^ 


1.46 X 10-^ 


+0.032 
-0.026 


'^■n'n+aa 


1.29 X 10"^ 


9.24 X 10-2 


+0.020 
-0.017 




2.53 


1.36 


+0.700 
-0.470 


-\-cc 


1.60 


8.61 X 10-^ 


+0.440 
-0.299 


^T]c-\-CC 


2.44 


1.31 


+0.680 
-0.450 




1.54 


8.29 X 10"^ 


+0.430 
-0.281 


O'ha + CC 


3.99 X 10"^ 


1.45 X 10"^ 


+0.165 
-0.089 


'^XaO + CC 


5.64 X 10"^ 


2.04 X 10-^ 


+0.232 
-0.128 


""Xel + CC 


6.25 X 10-^ 


2.26 X 10-^ 


+0.259 
-0.140 


''"Xe2 + CC 


2.44 X 10-^ 


8.90 X 10-2 


+0.100 
-0.055 



TABLE X: Total cross sections (in unit: pb) for the channel 
e+e" Z° IJ^cc) + ^ by setting = 1.50 ± 0.15 GeV. 
The final column presents the uncertainties from the center 
values under mc — 1.5 GeV. 



mi,(GeV) 


4.75GeV 


5.05GeV 


uncertainty 


0"T + 99 


6.95 X 10-2 


6.17 X 10-2 


+0.004 
-0.003 


0-T'+99 


3.44 X 10-2 


3.05 X 10-2 


+0.002 
-0.002 




8.72 X 10-2 


7.69 X 10-2 


+0.005 
-0.005 




4.32 X 10-2 


3.81 X 10-2 


+0.003 
-0.002 


""X + bS 


2.43 X 10-^ 


1.96 X 10-1 


+0.025 
-0.022 


<^r'+bb 


1.20 X 10-^ 


9.70 X 10-2 


+0.012 
-0.011 


"'Vh + bb 


2.17 X 10"^ 


1.74 X 10-1 


+0.023 
-0.020 


l^ri'^+bb 


1.07 X 10"^ 


8.61 X 10-2 


+0.011 
-0.010 


<^hb + bb 


8.12 X 10-=^ 


5.83 X 10"^ 


+0.001 
-0.001 


^Xb„ + bb 


1.28 X 10-2 


9.31 X 10-^ 


+0.002 
-0.002 


^XM +bb 


1.23 X 10-2 


8.87 X 10"^ 


+0.002 
-0.002 


^XHO+bb 


5.09 X 10-^ 


3.67 X 10"^ 


+0.001 
-0.001 



TABLE XI: Total cross sections (in unit: pb) for the channel 
e+e- Z° ^ + X by setting = 4.90 ± 0.15 GeV. 

The final column presents the uncertainties from the center 
values under mb — 4.90 GeV. 



uncertainties caused by varying tUc = 1.50 ± 0.15 
GeV are ~ 40% for the S'-wave case, which change 
to be ~ 70% for the P-wave case; when X = gg, 
the uncertainties for both the 5'-wave and P-wave 
cases are ^ 15% — 19%. 



• For the bottominum production, when X — bb, the 
uncertainties caused by varying nib — 4.90 ± 0.15 
GeV are ~ 11% for the S'-wave case, which change 
to be ~ 15% — 18% for the P-wave case; when 
X = gg, the uncertainties for both the S'-wave and 
P-wave cases are ^6%. 
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B. Properties of the color-octet processes 

1. Total cross sections 

As discussed in the introduction, we wih discuss the 
foUowing sizeable color-octet processes 

e+e-^Z" ^ \{QQ)[{l'S^,^'>),{l's[^'>)]g)+QQ 

ipQ + QQ , 

e+e-^Z" ^ |(Q0)[(l'^r)].9> +00^ /iQ + QO , 
e+e-^Z" ^ \{QQ)[{l's[^^)]g)+QQ^XQJ + QQ 
and 

e+e-^Z'' ^ \iQQ)[[l'S^,^^),[l's[^^)]g)+g , 
I 



and Table Km] shows 



And the ratios for the case oi X = g are similar. 

We make a discussion of the relative importance of 
different channels at the super Z factory with Ecm = ttiz 
and at the B factory with Ecm = 10.6 GeV : 

• At the super Z factory, the channel e+e~ Z'^ ^ 
\{QQ)[(l^ S[^^j]g) + QQ provides dominant con- 



e+e-^Z° ^ \{QQ)[{l'S^,^^)]g}+g^hQ+g , 
e+e-^Z" ^ \{QQ)[(l's[^^)]g) + g ^ XQJ + 9- 



It is noted that the total cross sections for the color-octet 
channels versus the collision energy have similar shapes 
with Figs. (|6|7|8|9p . 

We present their total cross sections at the e+e^ colli- 
sion energy Ecm in Tables IXU and IXIIII For the channels 
of X = cc or bb, Table |XlT] shows 



I 

tribution over other octet channels, due to the 
topology of Fig.(|3J;). While, at the B factory, the 
quarkonium production is dominant by the chan- 
nel via the 7* propagator, we will not have such 
big enhancement; i.e., we have 



a{e+e- ^ ^ \{cc)[{l\S^^^)]g) +X^J/ij + X) 



cr(e+e 



(T(e+e- ^ \ {cc)[{l^S\''>)]) + X ^ J/i' + X) 

aie+e- ^ Z° ^ \{cc)[{l'^s[^^)]g) +X^J/ij + X) 



ZO^\{ccm's[^^)])+X 



Z° 



\{bbm's^,'%)+x 



J/ij + X) 

^T + X) 



ZO ^ \{bb)[{l^si^'>)]) +X^T + X) 
^Z°^\{bb)[{l^s[^'>)]g)+X^T + X] 



cr(e+e 
(T(e"'"e" 



a{e+e- ^ Z" ^ \{bb)[{l^s[^'')]) +X^T + X 



= 0.1%, 
= 50%, 
0.01%, 
= 1.9% 



a(e+e- ^ Z° ^ |(cc)[(lig^^^)].g) +X^hc + X) 
aie+e- ^ ^ \{cc)[{l^ P^^^)]) + X ^ he + X) 
^(e+e- ^ Z" ^ \ice)[{l^S^^%) + X ^ XcJ + X) 
aie+e- ^ Z" ^ |(cc)[(l3pj^))]) + X ^ Xcj + X) 
aie+e- Z° ^ \ibb)[il' S^o^%) + X ^ h + X) 



aie+e- ^ Z^ ^ \ibb)[il^P^^'>)]) + X 



a(e^e ^■ 



Z''^\ibb)[il^s[^%)+X 



> hb - 
XbJ 



aie+e- Z^ ^ |(66)[(l3pj^')]> +X^ Xb.J + X) 



X) 
-X) 



2.9%, 
433%, 
= 5.5%, 
= 173%. 
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J/^ or T production channels 




cr(6) 


cr(c) 


O"tot 


e+e" Z° ^ Kcc)^:^^,^ +cc-> J/iP + cc 


1.65 X 10-3 


2.65 X 10-^ 




1.69 X 10-3 


e+e- ^ Z° |(cc)[(l='s{^')]g) + cc -> J/iJ: + cc 


1.71 X 10-3 


1.80 X 10-* 


8.92 X 10-* 


9.11 X 10-* 


e+e- ^ Z° ^ |(cc)[(l^S^'')]3) + 5 ^ J/^A + 








4.87 X 10-* 


e+e- ^ Z° ^ |(cc)[(l^Sr)]<?) + 5 ^ J/^A + .g 








3.31 X 10-3 


e+e" ^ |(66)[(l^S*^')]p) + 6&^ T + &6 


1.86 X 10-^ 


5.22 X lO-** 




2.76 X 10-^ 


e+e" ^ Z° -> l(6b)[(1^5*^')]5) + T + &6 


2.09 X 10-^ 


1.13 X lO-"^ 


4.05 X 10-3 


4.23 X 10-3 


e+e- ^ ^ + 5 ^ T + 5 








2.81 X 10-* 


e+e- -> Z« ^ l(fe6)[(l35P^)]p) +g^T + g 








5.93 X 10-* 



TABLE XII: Total cross sections (in unit: pb) for the color-octet heavy quarkonium production via the e+e annihilation at 
Ecm = Z*^ , where the subscripts (a), (6) and (c) stand for Feynman figures shown in Fie:s. (|3l4|l respectively. The symbol atot 

refers to the sum of cT(^a), <^{b); <^(c) and their interference terms for the channel e+e- — > — >■ | (QQ) [^1*Sq^''^ , |^l3S'j*^^]g) + 
QQ; and it refers to the cross section of Fig.© for the channel e+e- Z" ^ \ {QQ)[(l^ S^fA , (l^SfA]g) + g. 



hc{hb) or Xcjixbj) production channels 


cr(a) 


cr(b) 


O-(c) 


O"tot 


e+e- ^ ^ |(cc)[(l*S',',^')]g) + cc^ftc + cc 


6.59 X 10-3 


1.06 X 10-* 




6.78 X 10-3 


e+e- ^ |(cc)[(l3S'l^^)]3) + cc ^ xa + cc 


6.85 X 10-3 


7.21 X 10-* 


3.56 


3.65 


e+e- ^Z'^ |(cc)[(l*S<"')]5) + g ^ + g 








1.95 X 10-3 


e+e- ^Z°^ |(cc)[(l35l^')]g) + g ^ Xc.7 + <7 








1.32 X 10-^ 


e+e- ^ Z° |(&6)[(l*S',f + bb ^ ht + bb 


2.53 X 10-* 


7.09 X 10-^ 




3.75 X 10-* 


e+e- ^Z° ^ \{bb)[{l^S^^^)]g) + bb ^ Xbj + bb 


2.09 X 10-* 


1.53 X 10-* 


4.26 X 10-^ 


4.44 X 10-^ 


e+e- ^ -> j(b6)[(l*S<^')]3) + p ^ ht + g 








3.81 X 10-3 


e+e- ^Z°^ \{bb)[{l^S\^^)]g) + g ^ Xbj + g 








8.05 X 10-3 



TABLE XIII: Total cross sections (in unit: pb) for the color-octet heavy quarkonium production via the e+e annihilation at 
Ecm = Z*^ , where the subscripts (a), (b) and (c) stand for Feynman figures shown in Fie:s. (|3l4|l respectively. The symbol atot 

refers to the sum of o"(a), <^(b), c"{c) and their interference terms for the channel e+e- Z" \{QQ)[(^1^ Sq^^^ , (^^^ S[^^^]g) + 
QQ; and it refers to the cross section of Fig.© for the channel e+e- -> \{QQ)[(^1^ S^^^^^ , (^l^S[^^^]g) + g. Here the 

contributions to 3Po,3 3 from color octet component {QQ)[{1^ S[^'')]g) have been summed up. 



^(e+e - 



7* ^ \{cc)[{l^s[^^)]g) +cc^ J/^j + cc) 
>7* -> |(cc)[(l3S'(^))]) -f cc^ J/i/;-|-cc) 
7* ^ \{ccm's[^%) +cc^ xcj + cc) 



r ^ |(cc)[(l3py')]> + cc^xcj + cc) 

I 



3.1%, 



19%. 



• At the super Z factory, it is noted that the chan- 
nel e+e- ^ ^ \{QQ)[{1^ S^^^), (1^ Sf %) - 
is less important than the case of e+e — >■ Z° — ; 
|(QQ)[(1*S'^^'), {l^s[^'')]g)+QQ by at least an or- 



9 



der of magnitude for the production of tJjq, hg 
and xqj- While at the B factory, the channel 



7* ^ l(QQ)[l'^o^l9> + .9 is significant, 



a(e^ 



\(cc)[(l^sl,'%)+g^J/,p + g) 



a(e^ 



7* ^|(cc)[(l35p))])-f 



cc 



30%. 
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2. Dijferential cross sections for the octet channels 

To show the relative importance among different chan- 
nels as indicated by Tables IXIII and IXIIII more clearly, 
we present the Pt and rapidity y— distributions for the 
octet production channels of J/ip, T, hi,, XcJ and 
Xbj in Figs. (|15|16|) . Furthermore, the total cross sec- 
tions with different ptcut and ycut are presented in Ta- 
bles |XlV] and |XVl They show that the channel e+e^ 
Z° \{QQ)[il^s[^^)]g) + QQ will lead to a peak in low 
Pt region, about 5 GeV for charmonium and 10 GeV for 
bottomonium. 

As a final remark, to be useful reference, we present 
the charmonium and bottomonium Pt distribution with 
various rapidity cut [ycut) and their y distributions with 
various transverse momentum cut (Ptcut) for the chan- 
nel e+e- ^ Z° ^ \Hq) + QQ in Figs. (|l7ll8ll9l20l) . 
In these figures, both color-singlet and color-octet con- 
tributions are taken into consideration. Here, for the 
V'Q production through e+e" ^ ijjQ + QQ, the 

components \{QQ)[(\^ S^^^^]g) ,\{QQ)[{\'' sf^^]g) and 

l(QQ)[(^l'^'S'i"''''^]) have been summed up. While, for the 
P-wave quarkonium production through e+e^ — > Z'' — > 
{hQ,XQj) + QQ, the contributions to hq and XQJ Pro- 
duction, including both color-singlet and color-octet con- 
tributions, have been summed up. 



IV. SUMMARY 

In this paper, we have studied the charmonium and 
bottomonium production through the e+e^ annihilation 
in the leading ag-order. A detailed discussion on the 
heavy quarkonium production at the super Z factory 
via the two types of semi-exclusive channels, e+e^ — >■ 
\Hqq) + X with X = QQ or gg, have been presented. 

As usual, total cross sections are dominated by the 
color-singiet IS'-level quarkonium states. However the 
color-singiet quarkonium states at the 2S'-level and 
the IP-level, and the color-octet quarkonium state 

\{QQ)[^^ sf'^]g) can also provide sizable contributions. 
Total cross sections for the color-singlet channels ver- 
sus the collision energy Ecm have been presented in 
Figs. (|6l7l8l9p . The curves for the color-octet 5- wave 
states have the same shapes with the corresponding color- 
singlet iS-wave states. In the low collision energy region, 
the cross sections are dominated by the processes via the 
7* propagator (agree with observation at the B factory) : 
main contributions are around 6 — 20 GeV for charmo- 
nium production and 10 — 40 GeV for bottomonium pro- 
duction. Around the collision energy Ecm = friz, due 
to the Z°-boson resonance effect, the cross sections be- 
come much larger for the processes via the Z^ propaga- 
tor. Then, in even higher energies, the total cross sec- 
tion drops down logarithmically for all the processes. By 
varying Ecm within the range of (1 ± 3%)mz, the total 



cross sections changes down to of its peak value 

for e+e~ -> J/^p + cc, which is down to (Jg^) for 

e+e- ^ ZO ^ T + 66. 

At the super Z factory, the heavy quarkonium pro- 
duction processes e+e" — >■ \Hqq) + X arc dominated 
by those via the Z° propagator. Comparing to their 
productions at the B factories (running at the energy 
scale around 10.6 GeV), total cross sections of e+e" — )■ 
\Hqq) + X become much larger in the super Z-factory 
or the GigaZ program of ILC. A programming high lu- 
minosity of the super Z factory is also helpful for ac- 
cumulating data. For example, we estimate the num- 
ber of events under the luminosity C — 10'^'*cm~^s~"'^(~ 
10^pb~^/year [1^). Setting Nh to stand for the number 
of hadron H events generated per year, for the charmo- 
nium production, we obtain 

Nj/^ = 2.74 X 10^ ,N,^, = 1.16 x lO'^ ,7V^^ = 1.76 x 10^ 
iV,,- = 1.11 X 10^ ,Nh^ = 2.41 X 10^ ,N^^„ = 7.37 x lO", 
7Vx,i = 1.55 X 10^A^xe2 = 2.17 X 10^ 

for the process e+e^ |i?cc) + cc; 

Nj/^ = 3.91 X 10^ , N^, = 2.47 x 10^ , N,^^ ^ 1.72 x 10'', 
7V^^ = 1.09 X lO'' 

for the process e+e~ — )■ jiJce) +.9.9- For the bottomonium 
production, we obtain 

A^T = 2.22 X 10"* ,Nr' = 1.08 x 10"* , A^,,„ = 1.94 x lO'', 
N,-,'^ = 9.60 X 10^ , Nh, = 7.22 x 10^ , A^^.o = 1-58 x 10^ 

A^xw = 2.52 X 10^ , A^x,, = 2.90 x 10^ 

for the process e+e~ — >■ \H^£) + bb; 

Nr = 6.55 X 10^ , Nr' 3.24 x 10^ , 7V^, = 8.18 x 10^ 
= 4.05 X 10^ 

'lb 

for the process e+e~ — ?> l-ff^t) + 99- Thus, in addi- 
tion to the P-factories as BABAR and Belle and the 
hadronic colliders as Tevatron and LHC, the super Z- 
factory will provide an excellent platform for studying 
the heavy quarkonium properties. 

For the charmonium or bottomonium production chan- 
nels via the Z'^ propagator, there is approximate "spin 
degeneracy" . Taking the charmonium production chan- 
nels as an example, we obtain 

and 

'^e+e-^Z°^Xal+cc 

Such "spin degeneracy" is confirmed by analyzing the 
quarkonium pt- and y- distributions, and also by a cross 
check from the fragmentation approach. 
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FIG. 15: The charmonium (Left) and bottomoriium (Right) Pt-distributions via the production processes, 
|(QQ)[fl'S<,^') , fl^S$^')]3) + QQ^ {■^Q,hQ,XQj) + QQ- The colUsion energy E^^ = mz. 
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FIG. 16: The charmonium (Left) and bottomonium (Right) ^-distributions via the production processes, e^e — ^ 
\{QQ)\{\^ S^^^^ , (l^S'J^^^lff) + QQ (V'Q.'iQ.XQj) + QQ- The colUsion energy Ec^ = mz. 



We have also discussed the uncertainties caused by the 
effective quark masses. For the charmonium production 
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PtcutiGeV) 


5 


10 


15 


a{e+e~ ^ Z° \{cc)[{l'- S^o^^)]g) + a5 ^ J/ip + cc) 


1.62 X 10--^(96%) 


1.41 X 10-''*(83%) 


1.13 X 10-3(67%) 


a{e+e- Z° ^ \{cc)[{l^ S[^^)]g) + cc -> J/i; + cc) 


6.80 X 10^^(75%) 


4.41 X 10^^(48%) 


2.85 X 10^^(31%) 


a{e+e- Z° ^ \{cc)[{l''S^^^)]g) +cc^hc + cc) 


6.50 X 10"^ (96%) 


5.65 X 10"^(83%) 


4.54 X 10-^(67%) 


a{e+e- ^ Z° ^ \{cc)[{l^ Sf^)]g} + cc ^ XoJ + cc) 


2.72(75%) 


1.76(48%) 


1.14(31%) 


a{e+e- ^ \{bb)[(l'' S^o^^)]g} + bb ^ T + bb) 


2.69 X 10-^^(97%) 


2.46 X 10-^(89%) 


2.11 X 10-^^(76%) 


a(e+e- Z" \ {bb)[{i\s[^'')]g) + bh ^ T + hb) 


3.80 X 10"'' (90%.) 


2.90 X 10"-' (69%) 


2.02 X 10"-'(48%) 


a{c 1 c ^ Z" ^ |(W^)[(i'Sf + bb ^ -r bb) 


:.!.(.)(.) X 10 ^(97%:) 


X 10 '\mVC) 


2.87 X 10 -'(76%) 


a{e+e- ^ Z° ^ \{bb)[{V' Sf^)]g) +bh^ XbJ + bb) 


3.99 X 10-^(90%) 


3.05 X 10-^(69%) 


2.12 X 10-^(48%) 



TABLE XIV: Total cross sections (in unit: pb) for the octet-channels of charmonium and bottomonium production via 
propagator at the super Z factory with different Ptcut- mc=1.5 GeV and mi,=4.9 GeV. Here the percentages in the parentheses 
represent the ratios between the cross sections with various Ptcut and the one without Pt cut. 



Vcut 


0.5 


1.0 


1.5 


a{e+e- Z° ^ \{cc)[{l'- S^,^'')]g) + cc ^ J/iP + cc) 


6.45 X 10"'' (38%) 


1.18 X 10-3(70%) 


1.49 X 10-3(88%) 


a{e+e- Z° \(cc)[(l^ S^^^)]g) + cc ^ J /ij) -\- cc) 


4.26 X 10-1(47%) 


7.14 X 10-1(78%) 


8.48 X 10-1(93%) 


a{e+e- Z° ^ \{cc)[{l'- S^,^^)]g) + cc ^ he + cc) 


2.59 X 10-^(38%) 


4.73 X 10-3(70%) 


5.98 X 10-3(88%) 


a{e+e- ^ ^ \{cc)[{i' Sf%) + cc ^ XcJ + cc) 


1.71(47%) 


2.86(78%) 


3.40(93%) 


a{e+e- Z° l(66)[(li5'<^>)]g) + 66 ^ T + bfe) 


1.14 X 10-^(41%) 


2.07 X 10-^(75%) 


2.58 X 10-'5(93%) 


a{e+e- Z° ^ |(66)[(r'5f ')]g) + 66 -> T + 66) 


2.27 X 10-3(54%) 


3.63 X 10-3(86%) 


4.13 X 10-3(98%) 


a{e+e- Z° ^ \{bb)[{l^S^o^'')]g) +bb ^ ht + bb) 


1.55 X 10-'*(41%) 


2.81 X 10-''(75%) 


3.51 X 10-* (93%) 


a{e+e- ^ Z° ^ \{bb)[{l^S[^^)]g) +bb^ Xbj + bb) 


2.38 X 10-^(54%) 


3.81 X 10-^(86%) 


4.33 X 10-^(98%) 



TABLE XV: Total cross sections (in unit: pb) for the octet-channels of charmonium and bottomonium production via Z 
propagator at the super Z factory with different j/cut- mc=1.5 GeV and m6=4.9 GeV. Here the percentages in the parentheses 
represent the ratios between the cross sections with various ycut and the one without rapidity cut. 



Appendix A: Phase space splitting for the 
color-singlet case 

Generally, we can factorize the 2^3 phase phase into 
two parts 



d$3(Pl +P2;P3,P4,P5) 

(27r)V*) L+p^-j^Pi] n 

\ 1=3 / j=3 



{2n)^2Ej 



dg2 



= d^2iq;Pi,P2)^d^3{q;p3,p4,P5), (Al) 

where = (Ei=3^i)' " lEts^l' ft = {Ei,pi). 
At the same time, the hard scattering amplitude can be 
rewritten as 



-9,. + ] MU, 



A 



■^^^i-^i^3<(9,>)f-(9,A) 



= In M^iel{q,X^)\'\MUs<^4'l,^2)\\{A2) 

A1,A2 

where e{q, A) stands for the polarization vector. Then 

the process e~^e~ — \Hcc) + cc can be divided into two 
parts: the 2 — > 1 process (e+e" — > Z°) and the 1 — >■ 3 
process (Z° — )• \Hqq) + QQ). The phase space of 2 — >• 1 
process can be done easily. And the 1^3 process is the 
same as decays to three final particles. 
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